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We study the simultaneous determination of m s and V us from flavor-breaking 
hadronic r decay sum rules using weights designed to bring under better control 
problems associated with the slow convergence of the relevant D = 2 OPE series. 
Much improved stability and consistency is found as compared to the results of 
conventional analyses based on the "(fc, 0) spectral weights" . Results for m s are in 
excellent agreement with those of recent strange scalar and strange pseudoscalar sum 
rule analyses, as well as recent lattice analyses, while those for V us agree within errors 
with the output from recent lattice-based rpK^J/r^Tr^] and i^-based analyses. 
Very significant error reductions are shown to be expected, especially for V us , once the 
improved strange spectral data from the B factory experiments becomes available. 

PACS numbers: 12.15Ff,12.15.Hh,13.35.Dx,11.55.Hx 



I. BACKGROUND 

Measurements of inclusive flavor ij = ud, us vector (V) or axial vector (A) cur- 
rent induced hadronic r decay distributions yield kinematically weighted linear com- 
binations of the spectral functions, Pyh.^, of the spin J = and 1 parts, ^-yh.^, 
of the relevant current-current correlators. Explicitly, with Ry/A-ij = r[r~ — > 
v T hadrons y/j4;ii (7)]/T[t~ -> v T e~v e {i)\ Q, 

R v/A , l3 = 12n 2 \V l3 \ 2 S EW r ^ (l-y T ) 2 \{l + 2y T ) p^{s) - 2y T pf /A . ij {s 

Jth m r L 

where y T = s/m 2 ., V!y is the flavor ij CKM matrix element, Sew — 1-0201 ± 0.0003 jl] 
is a short- distance electroweak correction, and the superscript (0 + 1) denotes the sum 
of J = and J = 1 contributions. Since the spectral function combinations in Eq. ([1]) 
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correspond to correlators with no kinematic singularities, each term on the RHS can be 
rewritten using the basic finite energy sum rule (FESR) relation, 

dsw(s)p(s) = — - (p dsw(s)Il(s) . (2) 
2m J lsl=SQ 

Analogous FESR's, corresponding to spectral integrals, Ry/™:^, obtained by rescaling 

the kinematic weights in Ry/A;ij by (1 — y T ) k y™ before integration, are referred to as the 
"(fc, m) spectral weight sum rules". Similar spectral integrals and FESR's can be con- 
structed for So < m 2 , for general non-spectral weights w(s), and for either of the correlator 

combinations ^-v/A-ij(- s ) or s ^v /A -%j ( s ) • ^ e denote such spectral integrals generically by 
Rfj(so), and refer to the purely J = contribution in "inclusive" FESR's (those having 
both J = + 1 and J = contributions) as "longitudinal" , in what follows. 

V us and/or m s are extracted using flavor-breaking differences, 8R w (sq), defined by 

SR w (s ) = [Jft(ao) / \V ud \ 2 ] - [J£(*o) / |K S | 2 ] • (3) 

Since 8R w (s ) vanishes in the SU(3)p limit, its OPE representation, 5Rq PE (sq), begins 
with a dimension D = 2 contribution, proportional to m 2 . Experimental values for 
5R w (sq) over a range of so and w(s) allow m s and/or V us to be fitted, provided so is 
large enough that insufficiently-well-known higher D OPE contributions are small 0, 
Has a Si. As stressed in Refs. the smallness of m s (and hence of the 

flavor-breaking spectral integral differences) presents a challenge for the determination 
of m s but a significant advantage for the determination of V us . Explicitly, one has, from 
Eq. © 0, 

IK.I = ^RZ(^)/([R^o)/\V ud \ 2 } - 6R% PE (s )) . (4) 

At scales ~ 2 — 3 GeV 2 , and for weights used in the literature, the dominant D = 2 term 
in SRq PE (sq) is much smaller than the separate ud, us D = OPE contributions, and 
hence than the separate ud, us spectral integrals (for physical m s , typically at the few 
to several percent level). An uncertainty, A (5Rq PE (s )), in SRq PE (s ) thus produces a 
fractional uncertainty in \V US \, ~ A (5Rq PE (sq)) /2R™ d (so), much smaller than that on 
SRope( s o) itself. Moderate precision for 5Ro PE (so) thus suffices for high precision on 
\V US \, provided experimental errors can be brought under control. 

In what follows, we perform a combined extraction of m s and |V^ S | based on existing 
spectral data. The ud data [l(J [HI, 12] are already quite precise, but sizeable errors on 



the us data [13j, [14J, [15] limit currently achievable precision. We thus focus on better 
controlling uncertainties on the theoretical (OPE) side of the analysis, especially those 
associated with slower-than-previously-anticipated convergence of the relevant D = 2 
series 16 . 
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II. TECHNICAL ISSUES IN THE HADRONIC r DECAY APPROACH 

The first major stumbling block is the very bad behavior of the integrated longtitudinal 
D = 2 OPE series. Even at the maximum scale sq = m 2 allowed by kinematics, the 
series shows no sign of converging [l7j]. Even worse, all truncation schemes employed 
in the literature, with no exceptions, badly violate constraints associated with spectral 
positivity 0. A large part of the very strong unphysical /c-dependence seen in results 
for m s from inclusive (k, 0) spectral weight analyses was, in fact, associated with this 
spectral positivity violation problem Inclusive analyses employing the longitudinal 
OPE representation are thus untenable, and earlier results obtained using such analyses 
should be discarded in favor of those of more recent non-inclusive treatments 0, 0, E]- 

Fortunately, the severe problems of the longitudinal D = 2 OPE representation are 
easily handled phenomenologically, for a combination of chiral and kinematic reasons. 
Apart from the 7r and K pole contributions, longitudinal spectral contributions vanish 
in the SU(3)f limit and are doubly-chirally suppressed away from it. Since the values of 
the longitudinal kinematic weight at the K pole and in the region of the excited strange 
scalar and pseudoscalar (PS) resonances is essentially the same, this double chiral sup- 
pression is preserved in the ratio of non-pole to pole spectral integral contributions. The 
small residual non-pole us PS and scalar contributions can, moreover, be well-constrained 



phenomenologically, the former via a sum rule analysis of the PS channel |4|, [19J, the 
latter via i^7r-scattering-data-based dispersive analyses 0, 18, 20] (the most reliable be- 
ing the coupled- channel version discussed in Refs. |20j| , which incorporates short-distance 
QCD and chiral constraints). With the very accurately known n and K pole contribu- 
tions, these results allow a reliable bin-by-bin subtraction of longitudinal contributions 
to the experimental distribution, and hence a direct determination of the (0 + 1) spectral 
function, from which FESR's not afflicted by the longitudinal D = 2 OPE problem can 
be constructed. Note that agreement between the values of m s obtained from scalar and 
PS sum rule analyses employing the resulting us scalar and PS spectral "models" and 
those of recent Nf = 2 + 1 lattice simulations [22J rules out significantly larger non-pole 
longitudinal spectral strength. The residual non-pole part of the longitudinal subtraction 
is thus certainly small, and very well under control at the level required for (0 + 1) FESR 
analyses. We therefore focus, in what follows, on sum rules involving the flavor-breaking 
combination AIR» = 4°+^) - n£$™(a). 

The second problem concerns the slow convergence of the (0 + 1) D — 2 OPE series. 
For scales s ~ 2 — 3 GeV 2 , [An] OP£; is dominated by its D = 2 contribution 16], 



[AII(Q 2 )] ° D P ^ 2 = — ^ [1 + 2.333a + 19.933a 2 + 208.746a 3 + (2378 ± 200)a 4 + • • •] , 

(5) 

where a = a s (Q 2 )/n and m s = m s (Q 2 ), with a s (Q 2 ) and m s (Q 2 ) the running coupling 
and strange quark mass in the MS scheme. The 0(a 4 ) coefficient has been estimated 
using approaches previously successful in obtaining accurate predictions for the 0(a 3 ) co- 
efficient in Eq. ([5]) and n/-dependent 0(d 3 m 2 ) coefficients of the electromagnetic current 
correlator in advance of the explicit calculations of these values (23| . 
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Since, with 4-loop running, independent high-scale determinations of a s (Mz) [24j 
correspond to a(m%) ~ 0.10 — 0.11, Eq. §5§ shows that, at the spacelike point on \s\ = s , 
the convergence of the (0 + 1) D = 2 OPE series is marginal at best, even at the highest 
scales accessible in r decay. While \a s (Q 2 ) | decreases as one moves along the contour away 
from the spacelike point, allowing the convergence of the integrated series to be improved 
through judicious weight choices, this observation shows that, for weights not chosen 
specifically with this constraint in mind, one must expect to find very slow convergence 
of the integrated D = 2 series. The (k, 0) spectral weights, w^ k '°\y) = (1 +2y)(l — y) k+2 , 
withy = s/sq, are very much non-optimal in this regard, since |1— y\ = 2\sin((f)/2)\ (where 
<p is the angular position measured counterclockwise from the timelike point), is peaked 
precisely in the spacelike direction. Slow convergence, deteriorating with increasing k, is 
thus expected for the integrated D = 2 series of the (k, 0) spectral weights. The results 
of Table I of Ref. [l6| and row 1, Table I of Ref. [s| bear out this expectation. 

In evaluating the integrated (0 + 1) D = 2 OPE contribution, the size of the last term 
kept, the level of residual scale dependence, and the difference between the direct correla- 
tor and Adler function evaluations, both truncated at the same order, have all been used 
as measures of the truncation uncertainty. The slow convergence of the integrated series, 
however, can make it hard to be sufficiently conservative. For example, the quadrature 
sum of the last term size plus residual scale dependence, used previously in the literature 
to estimate the 0(d 3 ) Adler function truncation uncertainty, yields a result ~ 2.5 times 
smaller than the actual difference between the 0(a 3 )-truncated Adler function and (ob- 
truncated direct correlator results [1] . Since the growth of a s with decreasing scale makes 
higher order terms relatively more important at lower scales, premature truncation of a 
slowly-converging series typically shows up as an unphysical So-dependence in extracted, 
nominally s -independent quantities. With polynomial weights, w(y) = J2m c my m ^ f° r 
which integrated D = 2N + 2 OPE contributions not suppressed by additional powers of 
a s scale as c^/sq, such unphysical s -dependence can also result if unsuppressed higher 
D contributions which might in principle be present are incorrectly assumed negligible 
and omitted from the analysis. Since, typically, not even rough estimates of D > 6 
condensate combinations are available, such omission is most dangerous for w(y) having 
large values of the coefficients c m , with m > 2, where such D = 2m + 2 > 6 con- 
tributions will be enhanced. The (2,0), (3,0) and (4,0) J = + 1 spectral weights, 
w (2,o)( y ) = i^2y-2y 2 + % 3 - 7y 4 + 2y 5 , w^{y) = l-3y + 10y 3 - 15y A + 9y 5 - 2y 6 , 
and w^ 4 ' \y) = 1 — Ay + 2>y 2 + lOy 3 — 25y 4 + 24y 5 — lly 6 + 2y 7 , are examples of weights 
having such large, potentially dangerous higher order coefficients. 

In view of the above discussion, so-stability tests are essential components of any 
FESR determination of V us and/or m s . The existence of a stability window in sq for 
extracted quantities or, if not a stability window, then a window within which the ob- 
served instability is safely smaller than the estimated theoretical uncertainty, is crucial 
to establishing the reliability of the theoretical error estimate. 
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III. THE SPECTRAL AND NON-SPECTRAL WEIGHT ANALYSES 

The OPE and spectral integral inputs used in our analysis are outlined below. The 
fact that we restrict our attention to FESR's involving the V+A spectral combination, 
to weights satisfying w(s = sq) = 0, and to scales So > 2 GeV 2 , all serve to strongly 
suppress possible residual OPE breakdown effects [1, S, [I?]] . 

Integrals of the leading D = 2 contribution to [AH(Q 2 )] OPE , given in Eq. ([5]), are 
evaluated using two versions of the 0(a 4 )-truncated CTPT prescription [29j], one involving 
the RG-improved D = 2 correlator contribution [AIT]^^, the other the equivalent Adler 
function contribution, [AD(Q 2 )] ( ^ ) ^ = —Q 2 d[H(Q 2 )] c \^/dQ 2 . The difference between 
the two expressions, which would agree to all orders, but differ here by terms of 0(a 5 ) 
and higher, is used as one measure of the truncation uncertainty. For a and rh s we employ 
exact solutions corresponding to the 4- loop-truncated (3 and 7 functions [30], with initial 
condition a s (m 2 ) = 0.334 ±0.022. The remaining initial condition, m s (2 GeV), is either 
taken as input or obtained as part of the fit. 

Expressions for the D = 4 and D = 6 contributions are given in Ref . [l| . The dominant 
D = 4 term, proportional to (m s ss), is evaluated using ChPT quark mass ratios [28], 
GMOR for the light quark condensate, and r c = (mt££) / (m s ss) = 0.8 ± 0.2. D = 6 
contributions are estimated using the vacuum saturation approximation, and assigned 
an uncertainty of ±500%. Contributions with D > 6 are assumed negligible, s -stability 
studies being employed to test the self-consistency of this assumption. 

For the spectral integrals we work with the ALEPH ud and us [13j data, for which 
both data and covariance matrices are publicly available. A small global renormalization 
of the ud data is performed to reflect minor changes in the e, /i and total strange branching 



fractions since the original ALEPH publication. Following the prescri pti on of Ref. [31 
we also perform mode- by-mode rescalings to reflect current (PDG06 [25j ) values of the 
branching fractions of the various strange decay modes (32|. Errors on the K and 7r 
pole contributions are reduced by using the more precise expectations based on r[7r M2 ] 
and r[A" /t2 ]. With current data, ud and us spectral integral errors are at the ~ 0.5% 
and ~ 3 — 4% levels, respectively, for the weights to be discussed below. BABAR and 
BELLE will drastically reduce the size of the us errors in the near future. 



A. The (k, 0) spectral weight analyses 

As noted above, slow convergence of the integrated D = 2 OPE series is both expected, 
and observed jl, Ea], for the (k, 0) spectral weights. Further evidence that the OPE sides 
of the (k, 0) sum rules are not under good control is provided by Figure [TJ The figure 
shows the la contours for joint fits of m s and \V US \ to the so = m 2 experimental spectral 
integrals for various pairs of (k, 0) spectral weights. It is clear that no sensible common 
fit region exists, and hence that a reliable joint fit for m s and \V US \ cannot be obtained 
using the (k, 0) spectral weights considered previously in the literature. 

Further consideration is warranted for the (0, 0) analysis, which has been proposed 
in the literature as a particularly favorable one for the determination of \V US \ [t|. In- 
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deed, with m s from other sources as input, the So = m 2 . version of this analysis, would, 
if reliable, allow an improved determination of \V US \ with only improved us branching 
fractions as input, a feature special to this weight and this Sq value. Unfortunately, on 
the OPE side, the D = 2 truncation uncertainty appears much larger than previously 
anticipated [Sjj], as can be seen from (i) the rather poor agreement between the truncated 
correlator and Adler function evaluations; (ii) the strong deterioration in this situation 
with increasing truncation order; and (iii) the very significant so-instability in the out- 
put values of \V US \, at least for the PDG04 value of m s used as input in Ref. jsj. In 
Figure [2] we demonstrate that this so-instability is not an artifact of the particular m s 
employed in Ref. [j§]. The figure shows the OPE and spectral integrals as a function of 
so for the (0, 0) FESR, for a range of different fixed input m s . The value of \V US \ needed 
for this comparison was obtained by matching the OPE and spectral integral versions 
of 5R^°\m 2 T ). Note that very strong correlations exist amongst the OPE integrals for 
different So and, similarly, amongst the spectral integrals for different sq. Strongly dis- 
crepant so-dependences for the OPE and spectral integrals, as seen in the figure for the 
entire range of m s considered, thus rule out the possibility that an acceptable So-stability 
in \V US \ might be obtained from the (0,0) analysis, in its current form, for any input 
value of m s . Even restricting our attention to the rather narrow range of So values within 
0.4 GeV 2 of m 2 , the level of s -instability in \V US \ is > 0.0020, more than a factor of 2 
larger than previous estimates of the total theoretical uncertainty. In Ref. [8|, in order 
to obtain a more conservative estimate of the truncation uncertainty, we employed twice 
the sum, in quadrature, of the last included term and the integrated correlator-Adler 
function difference. The resulting estimated truncation uncertainty (total theoretical un- 
certainty) in \ V US \ was then ±0.0020 (±0.0022) for s = m 2 , compatible at least with the 
observed so-instability in the limited So region noted above. While this estimate is con- 
siderably more conservative than those employed previously in the literature, we argue 
that, given the observed level of s -instability, less cautious assessments can not be justi- 
fied. Since the D = 2 truncation uncertainty is unlikely to be reduced, we conclude that, 
unfortunately, theoretical uncertainties preclude a (0,0) spectral weight determination 
of \V US \ at the sub-±0.0020 level. Fortunately, alternate weights with significantly im- 
proved D = 2 convergence behavior exist, allowing one to take advantage of the general 
approach proposed in Refs. 0]. 



B. The non-spectral- weight analyses 

In Ref. [4j, three non-spectral polynomial weights, wio(y), wio(y), and W2o(y), designed 
specifically to produce improved integrated D = 2 convergence, were constructed. The 
significantly improved convergence is displayed explicitly in Ref. jij, and Table 1 of 
Ref. The weights were also designed to (i) keep higher order coefficients which 
might enhance D > 6 contributions small, and, (ii) strongly suppress spectral integral 
contributions from the region above s ~ 1 GeV 2 , where us errors are large. Here we 
consider also a fourth weight, w$(y), with lower degree, but less strong suppression of 
the high-s part of the us spectrum. While the latter feature leads to larger us spectral 
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integral errors with present data, the weight will be useful for future analyses based on 
data with smaller us spectral errors. Having an additional non-spectral weight also allows 
us to investigate more fully the issue of the mutual consistency of analyses associated 
with different improved-convergence non-spectral weight choices. 

In Figure [3] we display the matches between the OPE and spectral integral differences 
for W20 and Wio, for a range of input m s {2 GeV). The spectral integral differences are 
computed as described above for the analogous (0, 0) spectral weight case, shown in 
Figure [2J The W\q and w$ matches are intermediate in quality, and hence not shown 
separately. The contrast with the (0, 0) spectral weight case is immediately evident. 

In Figure @] we display the joint fit, s = m%., la contours, analogous to those of 
Figure HJ for various pairs of the non-spectral weights. Also shown are the la contours 
for two 3- fold and the full 4-fold fits. A good common fit region for m s and \ V US \ obviously 
exists, in sharp contrast to the situation for the (k, 0) spectral weights. 

We first consider results for \V US \ obtained, as in previous studies, using single weight 
analyses with sq = m 2 and external m s input. With the average of strange scalar and PS 
sum rule and Nf = 2 + 1 lattice results, m s (2 GeV) = 94 ± 6 MeV, advocated in the last 
of Refs. 0, one finds \V US \ = 0.2209 ± 0.0029 exp ± 0.0017 t/l for w 20 , 0.2210 ± 0.0030 exp ± 
0.0010^ for w 10 , 0.2206 ± 0.0032 e;cp ± 0.0007^ for w w , and 0.2218 ± 0.0037 ezp ± 0.0009^ 
for w$, with large experimental errors, dominated by those on the us distribution. The 
results are consistent, but cannot be directly averaged due to strong correlations. 

A cautious approach, given the instrinsically slow convergence of the (0 + 1) D = 2 
series, is to ignore external information on m s , perform a combined fit for m s and \V US \, 
and verify that the results for m s are consistent with what is known from other sources. 
The results for the 3-fold, s = w 2 o, Wi , w w fit, including theoretical errors on the 
same footing as experimental ones in the minimization process, are m s (2 GeV) = 89 ±26 
MeV and \V US \ = 0.2202 + 0.0046. The analogous full 4-fold non-spectral weight fit yields 
m s (2 GeV) = 96 + 31 MeV and |V US | = 0.2208 ± 0.0052. The larger errors in the second 
case result from the less strong suppression of the high-s us data region by w$. 

Since the combined joint fit results for m s are in excellent agreement with the external 
average noted above, it makes sense now to perform a combined so = m^. fit for \V US \ 
with m s {2 GeV) =94 + 6 MeV as input. The lower-error 3-fold fit (without w$) yields 

\V US \ = 0.2209 + 0.0031 , (6) 

which we take as our main result. Note that, had the fit been performed without including 
the uncertainty on the central value for m s , the output central value for \V US \ would 
have been shifted up to 0.2221. The result of Eq. (Q, though lower in central value, 
is compatible, within mutual errors, with both the recent determination, \V US \ = 
0.2249 + 0.0019 [33j (based on the conventional Leutwyler-Roos estimate for / + (0))) and 
the r[i^ 2 ]/r[7r M2 ] determination, \V US \ = 0.2223^o i 3 26 (based on the updated 2006 MILC 
evaluation of fx/fn |22|)- We comment further below on the status/reliability of the 
current hadronic-r-decay-based analysis. 
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IV. DISCUSSION, CONCLUSIONS AND PROSPECTS 

We have shown that it is possible to construct non-spectral weights which improve 
the convergence behavior of the integrated D = 2 OPE series and allow a self-consistent 
joint fit for m s and \V US \. While the results reported above for |V^ S | are compatible 
within errors with those from other sources, one should bear in mind that the errors 
on the present hadronic r decay data are sufficiently large that non-trivial shifts in the 
central value are likely when the B-factory data becomes available. To see why this is the 
case, note that the branching fractions for observed strange decay modes in the current 
data go down only to the ~ 3 x 10~ 4 level. Missing modes with lower branching fractions 
could have a non-trivial impact on the extracted value of \V US \. Missing a strange mode 
with branching fraction 1 x 10~ 4 , for example, would lower \V US \ by ~ 0.0004. At the 
desired level of precision, one thus needs to detect all strange modes with branching 
fractions at the few-to-several-10 -5 level. Such determinations are certainly feasible at 
the B-factory experiments where preliminary results for branching fractions at this level 
(with errors at the few-to-several-10 -6 level) have already been reported [34] . 

Regarding the goal of improved future precision, one should note that the ~ 1/2% 
uncertainty on the ud spectral integrals, though negligible in comparison to the us errors 
at present, contributes ~ 0.0005 to the uncertainty on \V US \, assuming the current as- 
sessment of the non-strange spectral errors is correct. One should also bear in mind here 
the disagreement between the isovector vector spectral contributions implied by isospin- 
breaking electroproduction data and those measured directly in r decay experiments |35| . 
For example, were the electroproduction results for -kit (whose implied r — > -k-kv t branch- 
ing fraction is 4.5<r below that measured directly in r decay) to be correct, the value of 
\V US \ obtained from the r decay analysis above would be raised by ~ 0.0018. Even the 
slightly lower central value of the preliminary BELLE r — > -ktiv t measurement (which is 
compatible, within its somewhat larger systematic error, with the earlier r world average 
value) would produce an increase of ~ 0.0006 in \ V US \. Resolving the discrepancy between 
the 7T7T electroproduction and r decay results is thus important not just for clarifying ex- 
pectations for the hadronic contribution to (g — 2) M in the Standard Model, but also for 
finalizing future r decay determinations of \V US \. 

We finally comment that, re-running the analysis with the central values of the us 
spectral distribution unchanged but with errors reduced by a factor of 5 (covariances by 
a factor of 25), and assuming m s (2 GeV) to be known to ±5 MeV, a combined total fit 
error below 0.0010 is achieved. Such an improvement in the us spectral distribution errors 
should be well within reach of the final B-factory analyses j36|. Further improvement on 
the theoretical component of the total error is almost certainly also possible through the 
use of new weights with significantly less strong suppression of the high-s region. Such 
weights, however, will become useful only once the errors in the us spectrum above the 
K* have been significantly reduced. The necessity of self-consistency and stability checks 
also means that an interaction between theorists and experimentalists will certainly be 
required for successful development of such weights. 
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FIG. 1: m s , \V US \ joint fits using sq 



and various pairs of (k, 0) spectral weights 
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FIG. 2: (0,0) spectral weight OPE and spectral integrals vs. so (GeV 2 ) for various input 
m s (2 GeV) and \V US \ obtained by matching the resulting OPE and spectral integrals at sq = m r 
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FIG. 3: u>20 and itfio OPE and spectral integrals (GeV 2 ) vs. so for various input m s (2 GeV) 
and \V US \ obtained by matching the resulting OPE and spectral integrals at sq = ntf. 



OPE vs Spectral Integral (w. 




OPE vs Spectral Integral (w 



Spectral Integral 
OPE 



(I.IKKi - 
0.004 - 
0.002 - 
(I - 



-0.002 I — 1 — ' — 1 — I — 1 — ' — 

1.5 2 2.5 3 3.5 1.5 2 2.5 3 3.5 



0.008 r 
0.006 - 
0.004 - 
0.002 - 
(I - 
-0.002 - 



1.5 2 2.5 3 3.5 1.5 



1.5 2 2.5 3 3.5 1.5 



2.5 3 3.5 





2.5 3 3.5 1.5 



2.5 3 3.5 



FIG. 4: 2-, 3- and 4-fold, sq = m 2 , joint m s , \V US \ non-spectral weight fits 
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